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The Monte Carlo model Jewel 1.0 (Jet Evolution With Energy Loss) simulates parton shower 
evolution in the presence of a dense QCD medium. In its current form medium interactions 
are modelled as elastic scattering based on perturbative matrix elements and a simple prescrip- 
tion for medium induced gluon radiation. The parton shower is interfaced with a hadronisation 
model. In the absence of medium effects Jewel is shown to reproduce jet measurements at LEP. 
The collisional energy loss is consistent with analytic calculations, but with Jewel we can go a 
step further and characterise also jet-induced modifications of the medium. Elastic and inelastic 
medium interactions are shown to lead to distinctive modifications of the jet fragmentation pat- 
tern, which should allow to experimentally distinguish between collisional and radiative energy 
loss mechanisms. In these proceedings the main Jewel results are summarised and a Monte 
Carlo algorithm is outlined that allows to include the Landau-Pomerantschuk-Migdal effect in 
probabilistic frameworks. 
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1. Introduction 

In ultra-relativistic heavy ion collisions, the produced QCD matter reduces significantly the 
energy of high transverse momentum partons. While most of the experimental evidence comes 
from studying the leading hadronic fragments of the parent partons via single inclusive hadron 
spectra and jet-like particle correlations at RHIC, there is ample motivation for studying medium- 
modified jets beyond their leading fragments, in particular: i) At the LHC a larger fraction of 
the entire medium-modified jet fragmentation pattern will become accessible above background, 
ii) Studying the distribution of subleading fragments is likely to discriminate between different 
microscopic mechanisms conjectured to underly jet quenching, thereby helping to characterise 
more precisely the properties of matter tested by jet quenching, iii) Modelling the distribution 
of subleading jet fragments is essential for any operational procedure aiming at disentangling jets 
from background or characterising the jet-induced modification of the background. Such reasons 
motivate the development of tools which account dynamically for the interaction between jet and 
medium, and which model medium-modified jets on the level of multi-particle final states. 

The Monte Carlo technique provides a powerful tool for the simulation of multi-particle final 
states. In the presence of medium effects, we expect that a parton shower can contribute to under- 
standing 'jet quenching', in particular for the following reasons: i) It reproduces the unmodified 
jet evolution as vacuum baseline, ii) Energy and momentum can be conserved exactly at each ver- 
tex, iii) Different microscopic mechanisms for the interaction between projectile and target can be 
tested, iv) The modelling of realistic multi-hadron final states is clearly beneficial for comparing 
theory and data. 

In developing Jewel (Jet Evolution With Energy Loss), our main focus was to arrive at a 
code which allows to study in detail the dynamics relating the evolution of the parton shower to the 
microscopic modelling of the medium interactions. 

This paper aims at summarising the most important features and results of Jewel 1.0. For 
more details and more elaborate discussion we refer to [[I], ||]. 

2. Vacuum Baseline 

In the absence of medium effects the Jewel parton shower is closely related to the mass- 
ordered shower in Pythia 6.4 [§]. It generates fragmentation patterns of a parton with energy E, 
produced in a hard scattering, into a multi-parton final state by probabilistic iteration of the 1 — > 2 
splitting processes. The virtuality Q 2 is used as ordering variable, since it is directly related to the 
parton lifetime. The probability that no splitting occurs between an initial and final virtuality Qf 
and Q 2 , respectively, is described by the Sudakov form factor 



Here, P a ->bc(z) are the standard LO parton splitting functions for quarks and gluons (a,b,c G {q,g}) 
and z is the energy fraction carried by the first daughter. The probability density for a splitting to 
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Figure 1: The thrust, thrust major and thrust minor (7J- = (T, T ma j, r m ; n )) distributions for ^/s = 200GeV 
e + e~ — > qq — > X collisions. Data of the ALEPH Collaboration [|]] are compared to simulations of Jewel. 
Dashed line: parton level after parton shower evolved down to Qq = 1 GeV; solid line: hadron level after 
parton shower evolution followed by hadronisation (Qq = 1 GeV). Figure taken from nil] . 



occur at Q 2 is given by 



s a (G?,e 2 ) 



d(ln(2 2 ) 



b,c 



where 



z + {Q\E) 



dz 



?-(e 2 ,£) 



q s (z(i-z)6 2 

271 



-A-*c(z)- 



(2.2) 



(2.3) 



One thus has to choose the parent parton's virtuality according to (|2.2j). The type of splitting 
is selected consistent with W a ^bc and the energy sharing is given by the splitting function P a ^bc(z)- 
This procedure is repeated for the daughters until no partons above the infrared cut-off scale Qq are 
left. Angular ordering is enforced by allowing only for splittings with decreasing emission angle. 

At <2o the parton shower is interfaced with a hadronisation model, that has to be flexible enough 
to be usable for the complex parton state in heavy-ion collisions. In Jewel a variant of the Lund 
string fragmentation [Q] is used, where the knowledge about the colour flow is replaced by the 
assumption of maximal colour correlation of partons close in momentum space. 

The Jewel parton shower is validated against a set of benchmark measurements in y/s = 
200 GeV e + e~ collisions at LEP. In Fig. [T] we compare the results of the simulation for the event 
shape variables thrust T, thrust major r ma j and thrust minor r m ; n to data, with the definitions, 
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(2.4) 



Thrust is a measure for how pencil-like an event is. For events, in which all momenta are (anti)parallel 
to the thrust axis Ht, T = 1 and for spherical events T = 1/2. Repeating the analysis in the plane 
transverse to the thrust axis one obtains thrust major. Thrust minor sums up the components of 
the final particle momenta p\, which are orthogonal to the plane defined by n and nj. 

As seen in Fig. [I], the final state parton shower provides a reasonable description of these jet 
event shapes over most of the measured range. The Jewel parton shower is not matched to exact 
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Figure 2: The inclusive distribution aW c h/d<i; with t, = In [Ejet/.Phadron] of charged hadrons in e + e~ — > — > 
X events at = 200 GeV. Data of the ALEPH Collaboration [|[] are compared to simulations of Jewel. 
Dashed line: parton level after parton shower evolved down to go = 1 GeV; solid line: hadron level after 
parton shower evolution to Qo = 1 GeV followed by hadronisation. Figure taken from [Q]. 



matrix elements, which may be the reason why it gives too few events with large l — T, T ma j and 
T m \a- Fig- % also shows that the thrust variables, since they are infrared and collinear safe quantities, 
are not very sensitive to hadronisation. 

The comparison to other event shape variables is of similar or better quality which shows 
that Jewel accounts with a sufficient accuracy for global features of the momentum flow. 

In comparison to jet event shapes, there are measurements which are more sensitive to the 
discrete and stochastic nature of the partonic processes underlying the QCD jet fragmentation. One 
such measurement is the n-jet rate, which is based on k± -clustering algorithms and quantifies the 
sub-structure of jets. The overall agreement of the simulation results with data is reasonably good, 
although the 3-jet rate is somewhat underestimated (possibly due to the missing matrix element 
matching). 

In contrast to the measurements discussed so far, the modelling of single inclusive intra-jet 
hadron distributions and multi-hadron correlations requires detailed knowledge about the hadroni- 
sation mechanism. This is seen for instance in Fig. [|, where we compare results of our simulation 
to data for the inclusive distribution (W c h/d£, with § = ln(£j et / ^hadron- Irrespective of the choice of 
<2o the partonic distribution is very hard while hadronisation is expected to soften the distribution 
considerably. As seen in Fig. ^[ Jewel with the modified string fragmentation describes the data 
to a level better than a few percent. 

We note that in the presence of a high-multiplicity environment, novel hadronisation mecha- 
nisms may play a role. Within these uncertainties, which are mainly related to the modelling of 
hadronisation, we have established that Jewel provides a reliable baseline for the characterisation 
of jet quenching phenomena. 

3. Medium Modification of Jets 

We regard the medium as a collection of partons acting as scattering centres. For the case of 
elastic interactions between the jet and the medium, each scattering centre displays to the partonic 
projectiles an elastic 2^2 scattering cross section da/df. For this exploratory study the medium 
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is modelled as an ideal gas of quarks and gluons with a constant temperature T. The masses of the 
scattering centres are fixed to m scatt = Hd(T) / V2, where Hd(T) is the thermal Debye mass. 

To specify the spatiotemporal structure of the parton shower, we start from the estimate that 
the parton shower evolves down to components of virtuality Qf on a time scale 1/2/. For a parton 
of energy E and mass Qf, this will be time-dilated in the rest frame of the medium to about E / Qj. 
If the parton originated from the branching of some parton of virtuality Q\, then the parton of 
virtuality Q t existed for a duration of approximately 



E E 

(3-1) 

/ 



For the parent parton, which initialised the parton shower, the lifetime is % = J^. In the case of a 
medium of constant density n scatterings occur during this time with the probability 

l-SnoscattO) = 1 -exp[-a e ias«Tj8] . (3.2) 

For the differential elastic scattering cross section, we choose the LO perturbative ^-channel parton- 
parton cross section, which gives the dominant contribution. To make contact with the previous 
studies, we have chosen two different regularisation schemes for the elastic scattering cross section, 
namely case I 

a- = / ' d\t\ ^I±M) Cr S ]+ {S :lf , (3-3) 



which is the default, and case II 
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Elastic scattering may occur between subsequent splitting processes and after the last splitting, 
provided the parton is still inside the medium. 

The model introduced here does not yet include a mechanism of radiative energy loss. Since 
this may be the main source of energy degradation, we have included an option to enhance the 
vacuum splitting functions by a factor (1 +/ me d), 

Pa^bc(z) — ► (1 +fmed)Pa^bc(z) , (3.5) 

as long as the splitting occurs in the medium. This prescription has been argued to display charac- 
teristics of radiative energy loss [Q]. 

It is assumed that the hadronisation of a jet is unaffected by the medium, i.e. occurs after all 
splittings and medium interactions. Since it is not clear, how jets in nuclear collisions hadronise, 
we investigate two options. One is that only the hadronisation of the parton shower is considered, 
whereas in the second the parton shower and the recoiling scattering centres hadronise together. 

In Fig. ||, we have calculated the resulting average energy loss for an in-medium path length of 
L = lfm in a medium of temperature T. Quantitatively, the models I and II show differences 
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Figure 3: The average parton energy loss dE/dx of a quark of energy E, undergoing multiple elastic colli- 
sions over a path length L = 1 fm in a thermal medium of temperature T. Elastic collisions are described by 
the infrared regulated partonic cross sections case I and case II. Figure taken from [[]]]. 
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Figure 4: Angle with respect to jet axis of recoiling scattering centres as compared to the undisturbed 
medium for different temperatures (cross section I, £j et = 100 GeV, L — 1 fm). Hadronisation is not included 
but may affect these distributions significantly. Figure taken from [[[[]. 



of approximately a factor 2 in dE/dx for a 10 GeV parton. Cross section II leads to a larger 
energy loss, as may be expected since there is minimum momentum transfer. We find that the 
temperature dependence of dE/dx shown in Fig. ^[ is consistent with the dependences reported 
previously |7]-[T3], p4| ). The differences between cases I and II are representative of the typical 
factor 2 uncertainties between different model studies, though some recent studies lead to slightly 
larger values of dE/dx than those shown in Fig. g[ see e.g. [13]. 

The energy lost by a jet is redistributed amongst the target components. Characterising the 
recoil of the medium may provide a means to disentangle different mechanisms of parton energy 
loss. In Fig. ^|the angular distribution of recoiling scattering centres is compared to the undisturbed 
thermal medium. The recoil moves predominantly in the jet direction with a maximum at A0 ~ 0.8 
nearly independent of the temperature. At face value, Fig. |] indicates that a jet can be accompanied 
by additional multiplicity which has its maximum separated from the jet axis by a characteristic 
finite angle A(j>. We note, however, that whether the partonic distributions of Fig. ^| will or will not 
change substantially upon hadronisation may depend on details of the hadronisation mechanism. 

Next, we study to what extent jet measurements are sensitive to medium effects. Only hadrons 
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Figure 5: Thrust, thrust major and thrust minor (T v = (T, r ma j, r m ; n )) for a single 100 GeV jet. The Jewel 
parton shower in the vacuum is compared to two scenarios including medium-induced parton energy loss. 
Left hand side: Collisional energy loss for a medium of T = 500 MeV and in-medium path length L = 5 fm 
(the recoil is hadronised with the medium). Right hand side: Radiative energy loss for / me( j = 3 and L = 5 fm. 
Only hadrons with energy above E cu i = 2 GeV are included. Figure taken from [|ljj. 



with energies above a soft background cut are included in the event shape analyses. This may indi- 
cate the extent to which different jet modifications remain visible above the hadronic background 
in heavy-ion collisions. Fig. || shows the thrust distributions for a single jet with and without 
medium modifications. Collisional energy loss is seen to lead to only a mild broadening of the 
distributions. The kinematics of elastic 2^2 scattering dictates that more energetic projectiles are 
deflected by smaller angles. The recoiling scattering centres have mostly relatively low momenta so 
that the medium-induced broadening becomes small if only hadrons of energy above Zs cut = 2 GeV 
are included in the analysis. In contrast to collisional energy loss, single partonic components of 
medium-induced additional radiation have a higher probability to carry a significant energy fraction 
of the initial projectile energy. So, on general grounds, one expects that the medium-induced broad- 
ening of the distributions in thrust, thrust major and thrust minor will persist even if soft hadrons of 
energy < E cut = 2 GeV are dropped from the analysis. This is clearly seen in Fig. § (right). The 
jet rates exhibit the same overall behaviour. Measurements of event shape variables and jet rates 
may thus provide a tool for disentangling collisional and radiative energy loss. 

In Fig. ||, we plot the medium modification of the single inclusive distribution dN/dE, in the 
presence of collisional and radiative medium effects. If elastic interactions with the medium are 
included but only the cascade is hadronised, then the total jet multiplicity of projectile hadrons 
depends only very weakly on the medium, since 2^2 processes do not increase the parton mul- 
tiplicity. However, the total jet multiplicity may increase significantly if recoil partons are counted 
towards the jet. Fig. |^ also indicates that radiative mechanisms result in a strong increase of intra-jet 
multiplicities. 



4. Outlook 



The most important future development of Jewel is to consistently include inelastic scatter- 
ing. In [16] a local and probabilistic implementation of the non-abelian Landau-Pomerantschuk- 
Migdal effect was presented. Analytic calculations [17] suggest that, in general, the LPM effect 
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Figure 6: The single inclusive distribution AN /dS, for a single medium-modified quark jet (E q = lOOGeV) 
after hadronisation (Qo = 1 GeV) (only charged hadrons are included). Collisional energy loss is calculated 
for T = 500 MeV and L — 5fm, with recoil partons either hadronised together with the cascade ('all') or not 
included in the hadronisation ('cascade'). Medium induced radiation is calculated for/ me( j = 3 andL = 5fm. 
Figure taken from [[!]] . 




Figure 7: The medium-induced energy loss AE of a quark of initial energy E q = 100 GeV as a function of 
the in-medium path length L. Results for the full MC algorithm with LPM effect and energy conservation 
(full case 1) are compared to: i) LPM effect without exact energy conservation (case 2), ii) incoherent limit 
without energy conservation (case 3), iii) incoherent limit with energy conservation (case 4), and iv) full 
case 1 with constraint that momentum transfer per scattering centre is limited to be soft (soft case 1). The 
bottom plot zooms into the small L-region. In this example L c ~ 9fm. Figure taken from [[To]]. 



can be implemented in a probabilistic Monte Carlo algorithm by requiring that the momentum 
transfer from different scattering centres to the partonic projectile acts totally coherently for gluon 
production, if it occurs within the formation time tf = 2co/k]_, and that it acts incoherently, if it 
occurs after tf. The Monte Carlo implementation of this prescription reproduces as a limiting case 
the BDMPS [18] results without modifying the gluon emission vertex. But it is more flexible in 
details of the modelling and has a natural way to implement exact energy-momentum conservation. 
In Fig. f7] the average energy loss for a projectile quark of initial energy E q = 100 GeV traversing 
a medium of length L is shown for different cases. Results for the full MC algorithm, including 
LPM-effect and exact energy-momentum conservation at each vertex are shown as 'full case 1'. 
Since the BDMPS multiple soft scattering approximation neglects the high-g^ tail of the scattering 
cross section, we also consider a case where the momentum transfer is limited to be soft ('soft case 
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1'). This case shows the characteristic L 2 -dependence of the average energy loss. Allowing for 
hard scattering, on the other hand, leads to a loss of coherence and thus increases the energy loss. 
'Case 2' is inspired by the high energy approximation assumed in many analytical studies. In this 
case the gluon energy has to be smaller than the initial quark energy, but the total radiated energy in 
unconstrained. As expected, this is a good approximation in the region L < L c ~ \/4a)^J7j, but for 
L > L c a quasi-incoherent regime is entered, with AE <=c L eventually violating energy conservation. 
'Case 3' is a scenario with incoherent gluon emission without energy conservation (but requiring 
(0 < E q ), whereas 'case 4' stands for incoherent gluon radiation with energy conservation. Energy 
conservation can be seen to become the dominant constraint for L> L c . The difference between 
cases 1 and 4 illustrates the quantitative importance of the LPM-effect. 

With this tool a dynamically consistent implementation of inelastic scattering including the 
LPM-interference in Jewel seems feasible. Among other planned future extensions are a realistic 
treatment of the medium and the generalisation to massive quarks. 
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